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Abstract 

The medium induced energy loss spectrum of a high energy quark or gluon traversing a hot 
QCD medium of finite volume is studied. We model the interaction by a simple picture of static 
scattering centres. The total induced energy loss is found to grow as L^, where L is the extent of 
the medium. The solution of the energy loss problem is reduced to the solution of a Schrodinger- 
like equation whose "potential" is given by the single-scattering cross section of the high energy 
parton in the medium. These results should be directly applicable to a quark-gluon plasma. 
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1 Introduction 



The determination of the radiative energy loss of a high energy charged particle as it passes 
through matter is a problem studied some time ago, in QED, by Landau, Pomeranchuk and 
Migdal There is recent data from SLAC |Q on radiative energy loss in QED [^]. New 

interest in this problem f^-^ has arisen because the corresponding problem in QCD, that of the 
energy loss of a high energy quark or gluon due to medium stimulated gluon radiation, may be 
important as a signal for quark-gluon plasma formation in high energy heavy ion collisions. 

Recently this problem was considered for infinite matter in P], hereinafter referred to as 
BDPS. Using the Gyulassy-Wang (GW) model for hot matter they observed that the QED 
and QCD problems are mathematically equivalent if one identifies the emission angle of radiated 
photons in the QED case with the transverse momentum of radiated gluons in the QCD case. 
The energy loss per unit length, —dE/dz, in hot QCD matter was found to be proportional to 
OsVE for an incident parton of energy E. The ^/E growth of —dE/dz was unexpected and 
suggested that energy losses of high energy jets in a quark-gluon plasma might be large. The 
procedure used by BDPS was not adequate for determining the exact logarithmic prefactor in 
—dE/dz. The problem was revisited in 0, hereinafter referred to as BDMPS, where a simple 
differential equation was given to determine the spectrum udl/dudz for radiated photons (in 
QED) and gluons (in QCD) which corrects the prefactor of the result in BDPS. 

The present paper generalizes the BDMPS approach to finite length hot matter. Although 
our discussion is carried out in the context of QCD it is a simple matter to change variables in the 
QCD results to get the corresponding QED results. Our main interest is in the situation where 
the incident quark or gluon is sufficiently energetic so that the length of the matter, L, satisfies 



L < Lcr = yXgE/jj,"^, with the gluon mean free path and /i the Debye screening mass of the 
medium. Our principal results are for udl /dudz, given in ( ^.16|) , and for dE/dz given in (|6.6|) 

and ( |6.7| ). The total energy loss AE in hot QCD matter of length L is AE = — — —^L^ ln(L/Ag) 

8 Xg 



as given in ( |6.8|) with R the colour representation of the incident parton. The total energy loss 
is found to be proportional to L^, a result which is surprising at first glance. However, by 
considering the limiting case where L = L^r it becomes clear that this corresponds to the BDPS 
result AE oc LcrVE. The fact that AE is proportional to has a simple physical interpretation. 
Since we expect the energy weighted spectrum ujdl /dujdz to be integrable in the infrared, AE 
is roughly determined by the maximum energy a radiated gluon can have still maintaining a 
coherence length < L. But, the formation time of a radiated gluon is 2ci;/A;^„^^, with k^rnax 
the maximum transverse momentum that the gluon gets by rescattering in the medium as it is 
being produced. Taking kj_^^^ = Lfi'^/Xg, with /i the typical momentum transfer to the gluon in 

a single scattering, and setting L oc . ^ , ^ , „ one finds u oc -^L^. This estimate, given more 

{L/Xg)lJ,^ 2 Xg 



precisely in (|6.9|), leads to the dependence of AE. 
The outline of the paper is as follows: 

In section ^ the emission probability of a soft gluon from a high energy quark traversing hot 
QCD matter is given in the GW model. The basic emission vertex is calculated in detail as 
are the subsequent rescatterings of the quark-gluon system passing through the matter as the 
gluon is becoming free. While, for simplicity, much of the discussion is given in the large-iVc 
limit, final formulas are given with exact colour factors. After observing that the QED and QCD 



2 



emission formulas are identical, with the identification of corresponding variables, the formula 
for the radiation intensity spectrum udl/dudz for infinite volume hot QCD matter is given as a 
direct consequence of the QED spectrum derived by BDMPS. 

In section | a heuristic discussion of energy loss is given both for L <^ Lcr and for L ^ Lcr- 
By requiring that these results match at Lcr a direct connection between AE oc L^J~E for L ^ L^r 
and IS.E oc l? for L <^ L^r is made. 

Section ^ is concerned with deriving the general equations governing radiative energy loss in 
a hot QCD plasma of extent L. The basic equation determining energy loss is a Schrodinger- 
like equation whose "potential" is given in terms of a single scattering cross section, in impact 
parameter space, of a high energy parton. We expect this same formalism, but with a different 
"potential", to apply to the energy loss problem in cold nuclear matter |Ty]. We presume, 
however, that the magnitude of the energy loss in hot and cold matter may be quite different. 

In section ^ we give an approximate solution for the radiation spectrum valid, at large L, 
for those gluon energies dominating the energy loss of the primary parton. Our main result is 
given in (5.16). We note that as L — oo the spectrum agrees with the result previously given in 
BDMPS. 

Formulas for the total energy loss due to medium induced radiation are given in section |^. 
We expect the most likely place that these results may have direct phenomenological application 
is in high energy jet production when a quark-gluon plasma is formed in heavy ion collisions. 
We also expect results much like (6.6) to hold in jet production in cold nuclear matter, but that 
is the subject for a subsequent paper |jTO|. 

Gluon emission with double scattering is given in Appendix A while rules for dealing with 
colour factors for multiple scattering are given in Appendix B. The calculation of the planar 
diagrams necessary to obtain (2.31) is given in Appendix C. A curious integral which arises in 
evaluating the energy loss is calculated in Appendix D. 

2 General expression of the medium induced radiation 
spectrum in QCD 

Here we derive the general form of the gluon energy spectrum induced by the propagation of a 
high energy quark in a finite non-abelian medium. Results are given for the case of an incident 
parton of arbitrary colour representation R. We also establish a formal analogy between QED 
and QCD. 

2.1 Model for multiple scattering 

In order to describe the successive interactions of a high energy incident parton with a hot QCD 
medium, we use the model introduced by Gyulassy and Wang and recently used in BDMPS to 
study the photon energy spectrum induced by multiple QED scattering of a fast charged particle. 
The main feature of the model consists in assuming that scattering centres are static. This allows 
one to focus on purely radiative processes, since the coUisional energy loss then vanishes. The 



3 



centre located at Xi creates a screened Coulomb potential 



H(f) = f , (2.1) 

47r \x — Xi\ 

with Fourier transform 

Vm = ^r^r—, e-^'-'' , (2.2) 

where g is the QCD coupling constant. We suppose that the range of the potentials Vi is small 
compared to the mean free path A of the incident parton, 

< A , (2.3) 

where /i is the Debye mass induced by the medium. This means that successive scatterings are 
independent, since the incident parton cannot scatter simultaneously off two distinct centres. 
As a consequence, its propagation is "time-ordered" , and we may number the scattering centres 
according to the interaction time (or equivalently the longitudinal coordinate) of the radiating 
parton. (See P] for a justification in terms of Feynman diagrams). Moreover, in the context of 
one gluon emission, this assumption allows us to neglect amplitudes involving four-gluon vertices 
of the type shown in Fig. 1. 



1 2 



Figure 1: This gluon emission amplitude induced by the scattering of an energetic quark off two 
static centres 1 and 2 is negligible in the A ^ limit. 

Finally, we work in the limit of very high energy E for the incident parton and in the soft gluon 
approximation, 

uj <^E . (2.4) 

Let us exhibit the implications of these two assumptions by giving the basic emission amplitude 
in a single scattering. 

2.2 Gluon emission induced by a single scattering 

The emission amplitude is depicted in Fig. 2. It includes the emission off the projectile (from 
now on chosen to be a quark) given by Mi and the emission off the exchanged gluon given by 
Ms. 

The colour indices of the static centre and of the incident quark are denoted by A, A' and B, B' 
respectively. The indices of the exchanged and radiated gluons are a and h. Neglecting screening 
for the moment, we write the amplitude for elastic scattering off a static source as 

M,, = n,sM%A ■ M%,A = g,,M''X'X ■ (2.5a) 
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Figure 2: Gluon emission amplitude induced by one scattering. 

Here 

Ml'^T = ^/(P. + P/)"^ '^0 ^I'A , (2.5b) 

where we neglected spin effects in the high energy hmit. The static source can be viewed as if it 
were a heavy quark. 

In Feynman gauge, the amphtude Mi (Fig. 2) for soft gluon emission may be expressed as 
the elastic scattering amplitude times a radiation factor as 

MiC^-gl^{T'T^)B'B-'j^{T^T')B'B]M%^ , (2.6) 
[k-pf k-pi J 

where e denotes the gluon polarization state. The generators of the fundamental representation 
of SU{N^) are T'^(a = 1, ... A^'^ - 1), satisfying [T", T^] = In the same way we get 

M2 - -g- ^ {g,,e ■ {pf - p,) - k,e, + k,e^] ■ M'^/X [T\ T^]b'b ■ (2.7) 

ypf - Pi) 

In addition to Mi and M2, there is a term M3 coming from gluon radiation off the static source. 
The sum of the three terms is gauge invariant. In a physical gauge such as light-cone gauge, M3 
is down by a factor of k^/uj compared to Mi and M2. In the calculation given below we use 
light-cone gauge and assume k^/uj <C 1. 



In a hot plasma the source is screened as indicated by (|2.1| ) and ( |2.2| ) in the GW model. The 
reader may have doubts as to the general gauge invariance of that model. These doubts may 
be put to rest by the following arguments. It is straightforward to show that Mi + M2 + M3 
remains gauge invariant when the emitted and exchanged gluons are given the same mass /x. As 
we shall see later, the emitted gluon has a small impact parameter for the physical problem we 
consider. As a consequence of the small impact parameter, one may neglect the mass for the 
emitted gluon; keeping the mass ji only for the exchanged gluon leads to the Gyulassy-Wang 
model. 

In light-cone gauge 

e = [eo, -eo, e^); e ■ k = ^ Eq = — — — ~ — . (2.8) 

00 + k// 2(jj 
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In the high energy hmit 



e ■ Pi e ■ pf ^ k± , , 

Li ^ lL ~ 2e± ■ — (2.9) 

k ■ Pi k ■ pf /c^ ' 



Ml ^ -2g el ■ ^ [T^T'^]B'B M^,^ . (2.10) 



where k± is the transverse momentum of the gluon with respect to the direction of the incident 
particle. Thus, 

kj 
kl 

In QED P|, the photon radiation amplitude vanishes in the limit E oo. In QCD, in the high 
energy limit only the purely non-abelian contribution to the gluon radiation spectrum survives. 
This is underlined by the presence of the commutator in ( |2.10|) . As a result we can use the 
eikonal approximation where the trajectory of the projectile is taken to be a straight line. Also, 

M, ^ 2g el ■ -^^ ~ f [^^ T%>s . (2.11) 

{k^ - qi_Y 

Finally, the radiation amplitude induced by one scattering of momentum transfer q reads 

Ml + M2 ~ -2g ■ J{k, q) [T^ T^^'b M%j, , (2.12) 
where the emission current J is defined as 

{k± - q±y 

We are interested in the gluon energy spectrum, which is given by the ratio between the radiation 
and elastic cross sections. Up to a common fiux factor 

- I 

da,, oc Cf\M-\^ dql ; Cp = , (2.14a) 

darad oc ^ CpN, P IM'^I^ d^kj^dq]_ , (2.14b) 

where jM^p is defined as iM^p^"'"^ = — Tr {WM"-'). Thus we obtain, for < < 

u:-^ = N,'^{j{k^qf) . (2.15) 
dud^k^ vr^ \ / 

As the amplitude has been evaluated for a fixed momentum transfer gl, an average over q± has 
to be performed. For this we use the probability density deduced from the elastic scattering cross 
section which is easily obtained from (PT^). Thus in ( [2.15] ) we define 



((...)) ^ Jd'q^V{qi){ ... ), (2.16) 
where the normalized cross section for elastic quark scattering reads 



7r(gi + /x2)2 
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with / c/^gl Vlq^) = 1. We have used the fact that the longitudinal transfer is negligible with 
respect to q± when E ^ oo. As we aim to derive the radiation density induced by multiple 
scattering, it is convenient to keep the colour structure together with the current and introduce 

Jeff = J{k,q)[T\T^]B'B . (2.18) 

The fact that the colour structure is the same as the three-gluon vertex allows one to give a 
compact diagrammatic representation of the effective current as shown in Fig. 3. 

k,b k,b 



Pi ,B ((^ pf, B' pi ,B ((^,Pf, B' Pi,B 



Pf 
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+ §1 ~^ k,b 

5 i fe 



Pi .B pf, B' 

k, b 



I 



I 

X 



a 



Figure 3: Diagrammatic representation of the effective current ( ^.18{ ). 
Then the differential energy spectrum is simply written as 



-^^ = ^V^ ■ (2-19) 

A comparison between ( |2.12| ) and ( p.l8| ) allows one to set the proper colour factor in order to 
normalize to the elastic scattering cross section. The square |</e//|^ includes the sum over all 
colour indices and the A^^^ in the denominator cancels the sum over initial quark colours while 
Cp corresponds to the colour factor of the normalizing elastic cross section. We see that the 
spectrum ( p^.l9D has exactly the same form as in QED IQ, up to the replacement of the photon 
angle by the gluon transverse momentum (and up to colour factors). 

The introduction of the effective current given in (|2.18D or in Fig. 3 will provide an important 



simplification in the case of multiple scattering. Let us indicate how this simplification appears 
in the case of two scatterings. 



2.3 Effective radiation amplitudes for double and multiple scattering 

Double scattering. For two scatterings, the radiation amplitude is given by a collection of 
seven diagrams. These are simply calculated in the framework of time-ordered perturbation 
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theory. We show in Appendix A that all amplitudes may be grouped in effective radiation 
amplitudes induced by momentum transfers gl or q2 at times ti and ^2 ; each is associated with 
a corresponding phase factor 



+J{k,q2)[c, h]a e'*^^ 

+J(A;-g2,gi)[[c,6],a] 6^*^^^+^*^ } (2.20) 



(see Appendix A for the notation concerning colour factors). This expression multiplies the 
elastic double scattering amplitude and may be represented diagrammatically as an effective 
emission current as in Fig. 3. This is shown in Fig. 4. 
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Figure 4: Radiation amplitude expressed in terms of effective currents. 



We thus have 

Mrad OC {Jeffik, q^)e'^^ + ljf{k, g2)e''^2 + Jeff{k - q^, qi) e'^^^ } , (2.21) 

where we use an obvious notation for the phases. The first term on the right-hand side of ( |2.21|) 
and Fig. 4 corresponds to gluon emission at ti followed by rescattering of the quark at t2- The 
second is gluon production at t2 while the third is gluon production at ti followed by rescattering 
of the gluon at ^2- As seen from ( p.21|) , quark rescattering does not affect the phase. 



Multiple scattering. The generalization of this simple result to N scatterings is straightfor- 
ward. After integrating over the time of emission t it is always possible to collect three pieces in 
order to construct the effective radiation amplitude induced by at time tj. Consider the three 
light-cone perturbation theory graphs of Fig. 5. 
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Figure 5: Diagrams contributing to the effective emission amplitude induced by the transfer qi±. 

For each diagram, integrating over t yields the difference of two exponential phase factors 
(see Appendix A). Keeping only the one depending on ti, we collect three terms having the same 
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phase, 



The sum of these three terms gives the effective current 



h I fc'^ 
> / "'i if L 

ki 
kl 

Ms - ^^^e^*4[c,a,] . (2.22) 



I k, c 

Iffik, q,) = J{k, q,)[c, a,] = •9J)J)J)J)Jb (2.23) 

/\ 

as in ( CT - 

Similarly to the case considered in (|2.2CI|) , the radiated gluon can rescatter on centres i + 
1,...N, so that the momentum k and the colour factor have to be changed accordingly in 
( p.23| ). For example, if the gluon emitted at centre i rescatters on centre j {j > i) the sum of 
the corresponding three terms results in an expression analogous to (|2.22| ), with k replaced by 
{k — Qj), since k labels the final real emitted gluon. In this case, one obtains 

Jeff{k - Qj, Qi) = J{k- Qj, Qi) aN-.-ttj+i ttj^i . . . ai+i[[c, aj],ai\ . (2.24a) 
This is diagrammatically shown below 



tj 

I I ! I 

mJ)]9J)MMJ)MJ)Jj)J)MS)S)S)JM^ 

J,jf{k-q„q,)= I I II I (2.24b) 

'dl \ \ \ \ Ij I 

X X X X X 

/ i+l j.i j N 

The associated phase is shifted according to 

exp{.4)^exp{.,(*^^.,^i^|-M} . ,2,5) 

In the total radiation amplitude, we should include, for centre i, the 2^~' possibilities (labelled 
by ir) for the quark-gluon system to rescatter on the remaining centres. For r = 1 to r = 2^^*, 
the associated phase gets modified each time the gluon rescatters (the phase is unchanged by 
quark rescattering) . Thus we write 

Mrad « E E q^)e''''^ , (2-26) 

1=1 r=l 

where the colour structure is included in Jeff- 
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2.4 Expression for the radiation spectrum induced by scatterings 

As for a single scattering, we square the radiation amplitude Mrad given in ( p. 261) and normalize 
by the multiple elastic scattering cross section to get the radiation spectrum induced by 
scatterings 

JJl!! = ^ f d'k, ( E f: %^ e^-^A . (2.27) 



du 7r2 J \^ ^ Nc Cp 



This expression deserves some comments. 




Figure 6: Interference term J^j-f ■ {Jiffy in the form of a connected graph. 



— 

1. In Jiff, the index i refers to the centre which induces the effective emission current. For a 
simple calculation of colour factors, it is convenient to represent interference terms in the 
form of connected diagrams, where the "conjugate amplitude" {JiffY appears in the lower 
part of the diagram (Fig. 6). 

2. The colour factor in the denominator of (|2.27|) corresponds to the normalization to the 
elastic scattering cross section, depicted in Fig. 7. 




Figure 7: Colour factor associated with the multiple elastic scattering cross section. 



This is easily calculated from the rules given in Appendix B. 

3. The sum over all possible gluon rescatterings is implicit in ( |2.27| ) (the sum over r in ( |2.2 
We should take into account all possible ways the quark-gluon system has to rescatter 
on centres n, in particular, for n > j + 1. However this part of the diagram describes 
the multiple scattering of the produced quark-gluon system after centre j, which has no 
influence on the energy radiation spectrum we are interested in. 
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4. Between centres i and j, it matters whether it is the quark or the gluon which absorbs the 
transverse momentum g^, because it changes the relative phase {ipi — (pj). To simphfy our 
derivation, we will first consider the large N^. limit, where all non-planar diagrams may be 
dropped, which corresponds to neglecting quark scatterings on centres n for i < n < j. 



5. Finally, the brackets in ( 2.27 ) denote the average over momentum transfers and longitudinal 
coordinates of the scattering centres. As in we have 

(( - )) ^ / n =^-P ("x) ■ / n A\ V{qlK ... ) , (2.28) 

where = ze+i — zi and V{q'j_) is the normalized cross section for elastic scattering, given 
by ( |2.17| ) in the case of Coulomb potentials. 



The products Jiff Jeff in ( p.27|) include a colour sum as indicated in ( p.30|) . We may rewrite 



as 



UJ- 



N N 



\ i=i j=i+l ^^c^F 



N 



As in the case of QED |^, we have the equivalent expression 



duj 



5 



N N j'i ji\ 

^eff ■ -Jeff (J(s>,-H>,) 



\ i=\ j=i+l 



e// ( 



1 + 



N ^ 
i=l 



(2.29a) 



(2.29b) 



In the large- iVc limit J\ff ■ J^jf is given by two sets of planar diagrams denoted by Y and H and 



shown in ( p.30[ ) 



Y 




mJ)J)J)^J)J)S)J)J)J)MJ)J)^M^id)j)^ 
II II 
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I 
I 
I 



i+l 




(2.30a) 



H 




I 
I 
I 

X 




(2.30b) 
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The second term of (p.29b| ) is the so-called factorization term, which corresponds to the limit 
of vanishing phases. In this limit, all emission amplitudes from the internal lines vanish (see 
Appendix A). Two cases have to be distinguished. 

• If the incident quark is produced at a time to = — oo, we see from Appendix A that only 
emission amplitudes from initial and final lines remain. The factorization term contribution 
is then equivalent to the contribution induced by a single scattering of momentum transfer 

(l±tot = H Qii ci-iid thus has a weak logarithmic medium dependence, as in the QED case 0]. 

i=l 

• In a realistic situation where the incident quark is produced, through a hard scattering, at a 
time to = 0, only emission from the final line remains (see the table of Appendix A). In this 
case the factorization term has no medium dependence at all, so that the medium induced 
spectrum is exactly given by the first term of ( 2.29b| ). It should be directly accessible 



experimentally by comparing hard scattering on a nucleus with that on a proton. 

We show in Appendix C that after dropping the medium-independent factorization term, 
( p.29b|) leads to the following medium-induced radiation spectrum in the large- iVc limit 



N N 



exp<jm^f//^ 



(2.31) 



i=l J=i+1 

where the t/j's are the transverse momenta of the gluon expressed in units of /i, 

Ui = Ui-i — Qi 

U=^;Q. = ^, (2.32) 



and Jj is the rescaled emission current 



J^ - . (2.33) 



The dimensionless parameter k is 



. = (2.34) 

2 u 

The radiation spectrum for the infinite medium QED case was given in for k ^ 1. We observe 
that the expression (|2.31| ) has the same form as in QED 0, with the replacements 



3 uj 1 

a — i> -Nc as : ^ — , (2.35a) 

photon "angle" Ui gluon "transverse momentum" Ui . (2.35b) 
This analogy allows us to give directly the result for the infinite QCD medium. Thus, 



.JL = 5^^Xl„i . (2.36) 
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Note that the radiation density is obtained by normahzing ( 2.31| ) to the distance A^A, where A 
is the mean free path of the incident quark ||^. 

The changes necessary to include all — corrections as well as the case of an arbitrary incident 

parton of colour represention R have been worked out in H]. These changes are 



A 



K 



2CfX 



2A 



9 ' 



(2.37a) 
(2.37b) 



This leads to the general formula 



UJ 



dl SagCn I, 1 
— ;- = ^ — y/«ln — 

diudz 27r An V k 



27r 



R 



\ 



XgU 



In 



UJ 



(2.38) 



where A^ is the gluon mean free path. 

We note that apart from the overall normalization, proportional to the squared colour charge 
Cr, there is no dependence of the induced radiation spectrum on the nature of the initial parton. 

We will now study, for a finite medium, the dependence of the radiation spectrum on the 
medium size L = NX. We present first a heuristic discussion of finite-size effects, and will use 
( p.31| ) as a starting point in section H. 



3 Heuristic discussion of the energy loss in finite length 
media 

When a very energetic parton of energy E is propagating through a medium of finite length 
L = NX the gluon radiation spectrum shows characteristic features depending on the gluon en- 
ergy UJ. For discussing the radiation density udl/du three different regimes may be distinguished 
PJ^ : the Bethe-Heitler (BH) regime with small gluon energies, the coherent regime (LPM) for 
intermediate uj, and the highest energy regime corresponding to the factorization limit. The 
coherent regime corresponds to the condition (cf. (4.19a) in |P) 

^<K<1 , (3.1) 



with K, given in ( p.34|) for the QCD case. Thus, for the finite media under consideration a 



reasonably large number of scatterings will be assumed, ^ 1. 

For the following qualitative derivations we neglect logarithmic factors. Thus we ignore 
numerical factors of order 1, and do not distinguish between propagating quarks and gluons. 
However we explicitly keep the parameters representing the medium. 



In terms of the gluon energy the condition (3.1) is 



UbH ~ A/i^ < < UJfact ~ — ^ < E . (3.2) 
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Obviously, ujfact < E only holds when L is less than the critical length, 



L< Lr 




(3.3) 



We note that the case uJfact <^ -E is consistent with the soft gluon approximation for the induced 
spectrum. 

The radiation spectrum per unit length behaves in the E ^ oo limit as 



UJ- 



dl 

dujdz 




L 



UJ < ujbh 

^BH < < iOfact 
^fact < UJ < E 



(3.4) 



for a finite length L < Lcr- These main features are illustrated schematically in Fig. 8. In the 
BH regime the radiation is due to N = L/X incoherent scatterings, whereas in the factorization 
regime the medium behaves as one single scattering centre. In the LPM regime l/y^ elementary 
centres act as a single scattering centre. 

In order to obtain the total energy loss AE we integrate the spectrum ( |3.4| ) over uj and z, 
with < UJ < E and < z < L. In addition to medium-independent contribution to the energy 
loss proportional to asE (a factorization contribution), we find the induced loss 



' A 



AE{L) ~ a, 



(3.5) 



It has been already pointed out in [11], that the total energy loss increases quadratically with the 
length L, and is independent of the parton energy E in the high energy limit. This interesting 
case is investigated in more detail in what follows. 

Here we conclude this heuristic discussion by considering the case of finite L, but with L > Lcr- 
This situation occurs for parton energies E < Ecr = L'^jj^/X (see (|3.3|)). By extending the 



coherent soft cu-spectrum 



up to ~ -E, the total induced loss is 



AE{L) 



Ifj^'^E 



X 



a, 



E Ec 



(3.6) 



which is i?-dependent and linear in L. The results given in ( p.5|) and (|3.6| ) are schematically 
summarized in Fig. 9, where the E dependence of AE is plotted for fixed L, and in Fig. 10, 
where E is fixed and the L-dependence is shown. 

The loss given by ( ^.6| ) is also relevant for an infinite medium, L — oo. An i?-dependent loss 
per unit length of propagating partons is found 



dE 
dz 



Ifi^E 



(3.7) 
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Figure 8: Radiation density l\3.4\) for a medium of finite length L < 
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Figure 9: Total induced energy loss as a Figure 10: Total induced energy loss as a 
function of the parton energy E. function of the medium size L. 



4 General equation for the induced radiation spectrum 

In this section we shall derive the general equations which govern the induced radiation spectrum 
for finite length materials. These equations generalize (4.34) and (4.40) of BDMPS. Our starting 
point, for the sake of simplicity, is the large-Nc formula (|2.31| ). We divide by L, the length of 
the material, and we allow the sum over scatterings to be arbitrary in number giving 



u dl ?>a^Nr 



dudz 2L7r^ 



/ rfA / d'^U ( 2Re^ Y'^Y.^i 

Jo J \ j=i+l { £=i 



r2^ 

X 



5 A - ^ A„ 



m=l 



(4.1) 

where A is the distance between the first and last scatterings. Expressing the sum over i as an 
integral of the position in the medium of the gluon emission vertex Jj and neglecting initial and 
final state scattering one arrives at 



UJ 



dujdz 2L7r2 
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2Re £ Ji ■ Jn+2 



n=0 



n+1 



A, 



(4,2) 



where now A is the distance between the emission vertices Ji and Jn+2- Although the sum in 
( [4. 21 ) formally goes from = to n = cxo we expect the typical number of scatterings to be L/\. 
It is convenient to scale all distances by the mean free path A. To that end let 



A_, 
A ' 



t — — , ti — , T — — . 



L 
A 



Then 



dl ?,a,N, 



I ^^'f^J n^"" ^ y n d'QeV{Q',)Ji ■ Jn+2 



oo „ n+2 



UJ 



dudz 2A7r Jo 

n+l 



n=0 " £=1 



r=l 



71+1 



exp < U^ts } - 1 



5 ft - E 

\ m=l / 



(4.3) 



where we have done the integral over Zi. In ([4 .31) the dependence on Qi and Qn+2 is contained 
only in the product of currents 



Ji ■ Jn+2 — -rpr 



Un+l — Qn+2 Un+l 



Ut iU,+Q,y UUn+l-Qn+2y U^+l 



(4.4) 



and in V{Q\) and ^(Q^+2)- As in BDMPS we integrate first over Qi and Qn+2 holding f/i, 
U2, ■ ■ ■ Un+l fixed. Defining an averaged elementary current 



foiUi) = I d'QiV{Ql)Ji =7i^ I dQ''V{Q 



we note that 

fo{Un+l) 

The spectrum can be written as 



d'^Qn+2 V{Qn+2) Jn+2 



dl 3a.Nr 



UJ- 



dudz 2tx\ 



2Re / dt{l- -) 



foiUi)-fiUi,t) 



K=0 



(4.5) 
(4.6) 

(4.7) 



where 



/(f/i,t) = /o(f/i)e 



00 n+l 



En/ d'QeViQ^^ 

n=l e=2 •' 

n+l / n.+l \ 

n / dU e-*^(^^^^^) t - ^ t„ f,{Un+i) . 

r=2 '' \ m=2 ) 



The first term on the right-hand side of ( [4.8|) is the j = i + l (nearest neighbours, n = 0) term 
of ( |2.31| ). The rest of the terms in ( [4.8| ) correspond to arbitrary numbers of gluon rescatterings 
between the currents 1 and ?t, + 2. 
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Equation (|4.8|) has a Bethe-Salpeter structure so it is straightforward to check that f{U,t) 
satisfies the natural integral equation 

/(?7, t) = f{U, 0)e-*(^-'"^') + J d^Q V{Q^) dt' e-(*-*')(i-^'^^')/(f7 - Q, t'), (4.9) 

where 

/([/,0)^ /;([/), (4.10) 

and where we have used 

Ui = Ut-i-Qi. (4.11) 

It is convenient to eliminate the t-integration on the right-hand side of ([4.9| ). This is achieved 
by taking a t-derivative which gives 

g-/iU,t) = -(1 -mf/2) fiU,t) + I d'Q V{Q') f{U-Q,t) ■ (4.12) 

Although ( [4.12| ) has been derived for an incident quark in the large- A^c limit, essentially the same 
equation holds for partons of colour representation R even without taking the large- A^^c limit. In 
the general case 

m, tR) = e-*«(i--^^) + ^ld^Q Vm f^dt'^ e-(*-*y(^— ^^)/([/-Q, 4) 

+ {^-^) r dt'a e-(*--*y(i--^^)/(f7, 4) (4.13) 



2CrJ Jo 

replaces (|4.9|) with = z/Xr, kr = ij?'\r/2uj and where \r = XCf/Cr. The second term on 



the right-hand side of (|4.13| ) corresponds to the emitted gluon rescattering which is dominant for 
incident quarks in the large A"c limit, and thus corresponds to the second term on the right-hand 
side of ( [4.9| ). The third term on the right-hand side of ( |4.13| ) is necessary to generate the non- 
planar graphs corresponding to the incident parton rescattering in the medium as discussed in 
BDPS. Indeed, as already discussed in section 0, there is no additional phase factor associated 
with the incident parton rescattering. As a result, this term is Q- independent and we can use 
Jd'QV{Q') = l. 

The colour factors appearing in ( [4. 13|) may be easily obtained from Appendix B. Taking a 
time derivative on both sides of ( 4.13| ), and defining 

A'c A'c . s 



we find 



«;=^«:, (4.15) 

—fiU, t) = -(1 - tKU')fiU, t)+ j d'Q ViQ') fiU -Q,t), (4.16) 

which equation is identical in form to ( [4.12| ). It is remarkable that the equation for the spectrum 
keeps the simple form 



(4.17) 
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provided we take with ( 4.16 ) the initial condition 

/(f7,o) = /o(t/) 



(4.18) 



instead of the initial condition implied in ( 4.13 ). 

We can make contact with the formalism of BDMPS for infinite length matter by taking 
To = cxD in ( [4.17| ). Identifying 

f{U)= / dTf{U,T) , (4.19) 



we see that ( 4.17| ) agrees with (4.28) of BDMPS when account is taken of the replacement 
3 

a — > -UgNc necessary to go from QED to QCD. Integrating ( [4.16| ) over r between r = and 



r = oo gives 



:i - ikU^) f{U) = /o(f7) + / V{Q^) f{U - Q) 



(4.20) 



which again reproduces the BDMPS result given in (4.30). 

Just as (4.20) was solved, for small k, by going to impact parameter space, it is convenient 
to recast ( [4.16| ) by defining 



which leads to 



d_ 



f{B,T)= -^J,Vl-{l-Vm) f{B,T 



(4.21) 



(4.22) 



(4.23) 



We can also use ( [4.21|) to rewrite the spectrum ( |4.171) as an integral over impact parameter with 
the result 



dl 3a.C 



TO 



2 Red dr 1 



dudz 27r^Ag | 7o V tqJ J 2tt B"^ 



d^B 1-^(^2 



-B-fiB,T] 



K=0' 



(4.24) 



the finite length generalization of (4.40) of BDMPS. In arriving at ( ^4.24| ) we have used 

J,{B) = -2m^^{l-V{B')). 



(4.25) 



Equations ( |4.23|) and ( [4.24| ), along with the initial condition / (-6,0) = / q{B) are the general 
equations governing the induced radiation spectrum for finite length matter in the context of the 
Gyulassy-Wang model of a hot QCD plasma. 



5 Solution for the radiation spectrum in the small k limit 



In this section we give an approximate solution to ( ^.23| ) and ( ^.24| ) valid in the small k limit. We 
restrict our evaluation of ( 4.24| ) to uj values not too far from those which dominate the energy 
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loss in finite length hot matter. The exact region of validity of our approach should become 
apparent as we proceed. 

Eq. ([4.23|) has a close resemblance to a Schrodinger, or diffusion, equation. However, since 
the general Schrodinger equation cannot be solved explicitely this analogy is not in itself useful. 
But, (|4.23|) is not only close to a Schrodinger equation it is even close to the Schrodinger equation 
for a two dimensional harmonic oscillator. To see this we recall that in infinite length matter the 
values of B which dominate the energy loss problem are those where B^ ~ ^/K, and we shall later 
see the same to be true in finite length matter. Also, recall that 1 — V{B'^) ^ jB"^ \n{l / B"^) for 
small B for Coulomb scattering. At small B^, ln(l/i?^) varies much more slowly than B^ so we 
may expect a reasonable approximation to consist in neglecting the variation of ln(l/i?^) when 
solving ( ^.231 ). In this case ( |4.23| ) can be solved in terms of the general solution to the Schrodinger 
equation for the 2-dimensional harmonic oscillator. In order to discuss the solution to (|4.23| ) in 
a context somewhat more general than that of Coulomb scattering write this equation as 



4/(5.r) 



f{B,T); f{B,0) = -^Bv{B') 



{5.1] 



where v{B'^) = 4(1 — V{B'^))/B'^. For potentials which decrease rapidly at large momentum 
transfer v{B'^) ^ v{0) at small B^ and v{0)fi'^ is just the mean momentum transfer squared. In 
the Coulomb case, v{B'^) ^ In 1/5^ for small B^ . We expect our approach to be valid whenever 



B' 



d 



dB' 



In v{B'^) < 1 



(5.2) 



for small B^. In the Coulomb case we expect our answer for udl / dudz to be correct within 
logarithmic accuracy in 1/k. 

So long as v{B'^) can be treated as a constant the solution to ( |4.23| ) proceeds in analogy with 
that of the two dimensional harmonic oscillator 



2m ^ 2 ° 



i,{B,t) 



(5.3) 



Comparing ( p.l| ) and ( p. 3D it is clear that one may take solutions to the harmonic oscillator if 
the identifications ^ 

m = — - , uo = V inv (5.4) 
2k 

are made. Thus we may immediately write the function / in terms of the Green function G as 

J{B, t)= [ d^B' GiB, r; B\ 0)7(5', 0) , (5.5) 



where 



GiB,r;B',0) 



mujQ 



27ii sin tuoT 



exp 



2 sin uqt 



B^ + B'^)cosuoT -2B ■ B' 



(5.6) 



with ujq and m as given in ( p.4| ). 
Eq. (^.5|) can be written as 



f{B,T) = N e-^^' / d^B' e 



-a{B'-Cf 



fiB',0) 



(5.7) 
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where 

a = — , J = -mujQtanuJoT , 6 = , N = ——. . (5. 

Using y/i = — i=r- one sees that Re a > and |a| ^ 1 for all values of r. Indeed for r <^ |c<Jo 



-1 



V2 



while for r ^ |c<Jo 



-1 



a ^ -— + —VT , (5.9a) 
Akt 12 



W— . (5.9b) 



4 

Hence |a|>l/v^ for all r. Thus so long as k is small the integral in (|5.7|) is dominated by values 
of \B' — (7p on the order of <C 1. Since we do not expect a rapid variation of / (B', 0) with 



B', the integral in ( ^.7| ) can be done by replacing B' by C in that function and then carrying out 
the gaussian integration over B'. One finds 

J(B,T) = —e-^^'JiC,0) . (5.10) 
a 

As we shall see below we always use values of r and u for which uqt is of order 1. Thus from 
(5.8) it is clear that C is of the same order as B. 
Recalling that 

f{C,0) = foiC) = --Cv{C') , (5.11) 



f{B,r) = - "^^(f exp I -'-muoB' tan uoA , (5.12) 
2 cos^ uJoT I 2 J 



we can write 

i7id{B^)B _ f i 

2 COS^ UJoT 

where we have used the fact that v is slowly varying to replace by B^ as the argument of 
that function. Finally, from the above it is apparent that B', in ( |5.7| ) is comparable to B. This 
occurs because the parameter k in (|5.1|) is small thus severely limiting the diffusion in B so that 

the time development of / {B, r) occurs with only moderate variations of B. That is, the impact 
parameter of the produced gluon is effectively frozen (at small values) during the whole time of 
its production. 

We can now use ( |5.12| ) in ( [4.241 ). The integration over r is made simple by noting that 
^ ^ 2iTiB d ( i o } 

J = ^exp|--mu;o5 tant^o^j • (5.13) 

We emphasize that in the case when {; is a constant, ( 5.13 ) is an exact solution to ( p.lj) . 
Substituting (|513|) into (|]24D leads to 

a;—— = - ^ Re H dr f d(BVB exp { --muoB^ tanuor} . (5.14) 
dujdz An^XgTo Jo J 12 J _ 
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We now do the S-integration, again neglecting the 5-dependence of v. One finds 



UJ 



dl 

duidz 



znXgTQ Jo 



tan uJoT 



(5.15) 



The 5-integration in ( |5.14|) is determined by values of on the order of |ma;o tanc^oTl"^. 
However, as is clear from ( p.l5D values of \ujQr\ ^ 1 do not contribute significantly to udl / dujdz, 
so we may assume that |ci;ot| > 1. This means that is on the order of Imc^ol"^ which, from 
(0) gives 52 . 



'k/v{B'^) as dominanting the spectrum of radiated gluons. In the Coulomb 
case this means B^ ~ ^Jk/ ln(l//t) while for potentials which decrease rapidly in transverse 
momentum, B'^ ~ is dominant. 



The time integral in ( |5.15|) is easily performed and gives 



u- 



dl GagCji 



dudz 



ttL 



In 



sm[uJoTo) 



(5.16) 



where we have set tqA = L the length of the matter. The coupling as in (|5.16|) should be 
evaluated at momentum scale on the order of fi'^/B'^ = fi^^Jiln 1/R)/R. We note that udl/doudz 
is small when \ujqTo\ <^ 1, while 



dl ?>asCR 

dudz |t>Joro|— >oo 2ti\„ 



(5.17) 



In the Coulomb case v{\/k) = |ln(l/fi:) bringing (|5.17|) in agreement with ( 2.38 ) and with (5.1) 
of BDMPS for infinite length matter. 



6 Energy loss in a finite length medium 

In order to find the total energy loss it is necessary to integrate ( ^.16|) over u. However, it is 
technically easier to go back to ( |5.15D and do the a;-integral before doing the r-integration. Thus 



dE 

dz 



dl SosCr 
uj- — —duj = — ^ Ke 



^0 dr r°° , 
— / duj2 



dudz 2TcXgTQ Jo T Jo 

It is convenient to change variables taking 

^ A/i^ _ 2x2 



UoT 



tan(c(JoT) 



Then 



K 



duj 



2uj 
1 



Ax 
I — 



Using ([6l^ ) in 



dE SagCR^"^ 

dz TTTo 



Re 



^ , 1°° dx 
vrdr / — 



'1 + i)x 



x^ [ tan[(l + i)x] 



(6.1) 



(6.2) 
(6.3) 

(6.4) 
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Using (see Appendix D) 



we finally arrive at 



Re 



X- 



dx 

~3 



{l+i)x 
tan[(l + i)x] 



TT 

6 



(6.5) 



— = ^^/^ ^0 v{1/to) = -jL v{1/to) . (6.6) 



In determining the argument of v we note that x is of order 1 in ( |6.4|) while r is of order tq. From 
( |6.2| ) we have that k is of order [tq?;]~^ in the dominant region of integration giving of order 
['i^o'^]"^ the effective argument of v. Thus in the Coulomb case and with logarithmic accuracy 

dZ 8 \g Xg 

The total energy loss in the matter is just L times ( p.7|) giving 

-AE = ^^LMn^. (6.8) 

8 Xg Xg 

We note here, as was noted earlier by BDPS, that the energy loss depends on the hot matter 
only through the parameter /i^/A in the logarithmic approximation. The fact that the total 
energy loss scales as is remarkable and, to our knowledge, has not been anticipated by earlier 
discussions in the literature. This loss in energy is due to emission of gluons having typical energy 

u = {XrrC^^^^^^^^^ . (6.9) 

^ ' X 4:X^ Xg 160 ^ ' 

In arriving at ( |6.9|) we have taken r = as the median value of r in (|6.4] ) and we have set 

X ~ 3 (see p.5|) as the median value of x in (|6.5|) . 

Of course any estimate of energy loss in possibly realistic circumstances in heavy ion colli- 
sions is hazardous and should be received with caution and scepticism. Nevertheless, we feel it 
necessary to say something about the size of the results (|6.7|) and (|6.8|). For hot QCD matter 



having temperature T = 250 MeV, an estimate |TJ], using simple perturbative formulas, gives 
A ^ Ifm and fi^/X ^ IGeV/fm^ For L = lOfm we find from (§]§ that AE = 80a, GeV while 
for L = 5fm, AE = 20as GeV, which are rather large numbers if as>l/3-l/2 or so. Note also 
that, from ( |6.9| ), typical values of u for L = lOfm are of order 3 GeV, which means that several 
soft gluons are emitted in a typical event. 
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A Radiation amplitude induced by two scatterings 



As a means of producing an energetic quark we consider deep inelastic scattering where the 
quark is produced at a time to- H = ~oo we have the circumstance considered by BDPS and 
BDMPS, while to = corresponds to a quark produced in the medium. 

The radiation amplitude is given by the sum of the seven diagrams shown below, with ti and 
t2 the interaction times of the incident quark or of the radiated gluon with static centres 1 and 2. 
The emission time of the gluon is denoted by t. The radiation amplitude is simply calculated in 
the framework of time-ordered perturbation theory. For each diagram, we indicate the associated 
phase factor and the interval over which the emission time t has to be integrated. The phase 
is given by the energy difference between the states after and before the interaction. We then 
exhibit the resulting difference of phase factors, the corresponding kinematical factors involving 
transverse momenta (see (|2.1CI|) and (|2.11|) ) obtained as in the case of a single scattering, together 
with the associated colour structure. A colour factor (T^T°'T'^) b' b will be denoted as bac. We 
also factor out the elastic double scattering cross section. 




to h t2_ 



I I 
I I 
I I 

X X 



oc / dt e 2- 
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(A.l) 
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7^ bca 



(A.2) 
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X X 



fc2 )2 / (fe-q2)i ('=-92)1 \ f/tl —09)1 

When we take the production time to be = —00, exponential factors containing to disappear 
and the radiation amphtude induced by two scatterings takes the simple form given in (|2.20|) . 

B Diagrammatic calculation of colour factors 

The generators T° (a = 1, . . . A''^ — 1) of the fundamental representation of SU{Nc) satisfy 

[T^T''] = ip^'^T" , (B.la) 

and 

[t\T^] = 1-5''^ + (f^^ , (B.lb) 

where the structure constant f"-^^ and the symbol d"-^'^ are respectively totally antisymmetric and 

symmetric. 

Using 

1 

(B.2a) 

ik , (B.2b) 

(B.2c) 

together with the decomposition 

T^T^ = 1 [T", T^] + ^ {t^ T^] , (B.3) 

and the relations [p!^ , p!^] 

rpbrparpb ^ ^rpa ^ X 

2iVe ^ ' 

jaib jhjc jcka _ _^^j:ijk ^j^^ 

2 



TriT^T^) 


= Tn 


5"^ = 


2^ ' 








2iV, ' 


= Cf 


^ik = 


jacd jbcd 


= Ca 


5'^* = 
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we get the following simple diagrammatic rules: 



nrsTflinnrc 




(B.5) 



(B.6) 



(B.7) 



1 



2 



1 

i 



(B.8) 



1 



I 

i 



(B.9) 



2 



(B.IO) 



C Expression of the medium induced spectrum 



We start from the first term of (|2.29b|) and calculate the product of currents given there and 



in ( p.30 ). After averaging over the trajectory of the incident quark, the contribution of the 
interference Ji ■ Jj e*'-'^'~'^j^ will depend only on (j — i). Thus we simply evaluate this product 
for i = 1 and j = n + 2, where n = j— i — 1 is the number of intermediate gluon rescatterings. 
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Since the probability density V{q^) is isotropic, we can choose the momentum convention 
separately for graphs Y and H as shown in ( p.30| ). We separate the colour and momentum 
structure of these graphs as 

Y + H = YcYm + HcHm . (C.l) 

We get 

Ym = /i ^Jl ■ Jn+2 = —Hm , (C.2) 

where the dimensionless current Jj is given in ( |2.33| ) in terms of the dimensionless momenta 
( p.32| ), and, from Appendix B, 



3l 

2Cp 



n+l 



-2H, 



c 



(C.3) 



where the factor NcCp'^'^ in the denominator accounts for the normalization to the multiple 
elastic scattering cross section. 

Finally the phases associated with Y and H are equal. 



/ tt2 n+2 Tj2 _ Tj2 ' 

^ - t£ 2 ITT\ 



n+2" 



Un+2 
2UJ 



A 



(C.4) 



Collecting (|C.1|) to (|C.4|) we get the large- iV^ formula (p.31|). 



D A curious integral 

Here we evaluate the integral (|6.5|) , 



g{x) dx 



Re- 



'1 + i)x 



tan[(l + i)x] 



1 



sh 2x + sin 2x 

X— 1 

ch 2x — cos 2x 



(D.la) 
(D.lb) 



by using the integration contour in the complex z = x + iy plane as shown in Fig. 11. 
Adding the contributions from the paths Ci and C3 gives for e 



Ici +Ic3 = '^ lini / g(x) dx = 21 



(D.2) 



When closing the path at infinity the contribution along C2 vanishes, whereas due to the single 
pole of the integrand at z = we have 



27ri 



'C4 



TT 

3 



:d.3) 



Observing that the poles of 1/ tan[(H-2)z], Zn = |(1 — i)nn are not in this quadrant we find the 
result given in (|6.5| ) for J, 

. \ . . 1 . TT 



'C4 



6 



:d.4) 



27 



3^ = Imz 




X = Rez 



Figure 11: Integration contour for the integral 



Finally, we evaluate numerically the median value of x in (p. la]) defined by 



g{x) dx = - 1 = — . 
2 12 



We find 



(D.5a) 



•^med — 3.2 



(D.5b) 
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